Overview of the Americas' First Peopling from a Patrilineal Perspective: New Evidence from the Southern Continent.

Uniparental genetic systems are unique sex indicators and complement the study of autosomal diversity by providing landmarks of human migrations that repeatedly shaped the structure of extant populations. Our knowledge of the variation of the male-specific region of the Y chromosome in Native Americans is still rather scarce and scattered, but by merging sequence information from modern and ancient individuals, we here provide a comprehensive and updated phylogeny of the distinctive Native American branches of haplogroups C and Q. Our analyses confirm C-MPB373, C-P39, Q-Z780, Q-M848, and Q-Y4276 as the main founding haplogroups and identify traces of unsuccessful (pre-Q-F1096) or extinct (C-L1373*, Q-YP4010*) Y-chromosome lineages, indicating that haplogroup diversity of the founder populations that first entered the Americas was greater than that observed in the Indigenous component of modern populations. In addition, through a diachronic and phylogeographic dissection of newly identified Q-M848 branches, we provide the first Y-chromosome insights into the early peopling of the South American hinterland (Q-BY104773 and Q-BY15730) and on overlying inland migrations (Q-BY139813).

Weaving Mitochondrial DNA and Y-Chromosome Variation in the Panamanian Genetic Canvas.

The Isthmus of Panama was a crossroads between North and South America during the continent's first peopling (and subsequent movements) also playing a pivotal role during European colonization and the African slave trade. Previous analyses of uniparental systems revealed significant sex biases in the genetic history of Panamanians, as testified by the high proportions of Indigenous and sub-Saharan mitochondrial DNAs (mtDNAs) and by the prevalence of Western European/northern African Y chromosomes. Those studies were conducted on the general population without considering any self-reported ethnic affiliations. Here, we compared the mtDNA and Y-chromosome lineages of a new sample collection from 431 individuals (301 males and 130 females) belonging to either the general population, mixed groups, or one of five Indigenous groups currently living in Panama. We found different proportions of paternal and maternal lineages in the Indigenous groups testifying to pre-contact demographic events and genetic inputs (some dated to Pleistocene times) that created genetic structure. Then, while the local mitochondrial gene pool was marginally involved in post-contact admixtures, the Indigenous Y chromosomes were differentially replaced, mostly by lineages of western Eurasian origin. Finally, our new estimates of the sub-Saharan contribution, on a more accurately defined general population, reduce an apparent divergence between genetic and historical data.

Archaeogenomic distinctiveness of the Isthmo-Colombian area.

The recently enriched genomic history of Indigenous groups in the Americas is still meager concerning continental Central America. Here, we report ten pre-Hispanic (plus two early colonial) genomes and 84 genome-wide profiles from seven groups presently living in Panama. Our analyses reveal that pre-Hispanic demographic events contributed to the extensive genetic structure currently seen in the area, which is also characterized by a distinctive Isthmo-Colombian Indigenous component. This component drives these populations on a specific variability axis and derives from the local admixture of different ancestries of northern North American origin(s). Two of these ancestries were differentially associated to Pleistocene Indigenous groups that also moved into South America, leaving heterogenous genetic footprints. An additional Pleistocene ancestry was brought by a still unsampled population of the Isthmus (UPopI) that remained restricted to the Isthmian area, expanded locally during the early Holocene, and left genomic traces up to the present day.

Y-chromosome and Surname Analyses for Reconstructing Past Population Structures: The Sardinian Population as a Test Case.

Many anthropological, linguistic, genetic and genomic analyses have been carried out to evaluate the potential impact that evolutionary forces had in shaping the present-day Sardinian gene pool, the main outlier in the genetic landscape of Europe. However, due to the homogenizing effect of internal movements, which have intensified over the past fifty years, only partial information has been obtained about the main demographic events. To overcome this limitation, we analyzed the male-specific region of the Y chromosome in three population samples obtained by reallocating a large number of Sardinian subjects to the place of origin of their monophyletic surnames, which are paternally transmitted through generations in most of the populations, much like the Y chromosome. Three Y-chromosome founding lineages, G2-L91, I2-M26 and R1b-V88, were identified as strongly contributing to the definition of the outlying position of Sardinians in the European genetic context and marking a significant differentiation within the island. The present distribution of these lineages does not always mirror that detected in ancient DNAs. Our results show that the analysis of the Y-chromosome gene pool coupled with a sampling method based on the origin of the family name, is an efficient approach to unravelling past heterogeneity, often hidden by recent movements, in the gene pool of modern populations. Furthermore, the reconstruction and comparison of past genetic isolates represent a starting point to better assess the genetic information deriving from the increasing number of available ancient DNA samples.

Analysis of the human Y-chromosome haplogroup Q characterizes ancient population movements in Eurasia and the Americas.

BACKGROUND: Recent genome studies of modern and ancient samples have proposed that Native Americans derive from a subset of the Eurasian gene pool carried to America by an ancestral Beringian population, from which two well-differentiated components originated and subsequently mixed in different proportion during their spread in the Americas. To assess the timing, places of origin and extent of admixture between these components, we performed an analysis of the Y-chromosome haplogroup Q, which is the only Pan-American haplogroup and accounts for virtually all Native American Y chromosomes in Mesoamerica and South America. RESULTS: Our analyses of 1.5 Mb of 152 Y chromosomes, 34 re-sequenced in this work, support a "coastal and inland routes scenario" for the first entrance of modern humans in North America. We show a major phase of male population growth in the Americas after 15 thousand years ago (kya), followed by a period of constant population size from 8 to 3 kya, after which a secondary sign of growth was registered. The estimated dates of the first expansion in Mesoamerica and the Isthmo-Colombian Area, mainly revealed by haplogroup Q-Z780, suggest an entrance in South America prior to 15 kya. During the global constant population size phase, local South American hints of growth were registered by different Q-M848 sub-clades. These expansion events, which started during the Holocene with the improvement of climatic conditions, can be ascribed to multiple cultural changes rather than a steady population growth and a single cohesive culture diffusion as it occurred in Europe. CONCLUSIONS: We established and dated a detailed haplogroup Q phylogeny that provides new insights into the geographic distribution of its Eurasian and American branches in modern and ancient samples.

Reconstructing the genetic history of Italians: new insights from a male (Y-chromosome) perspective.

BACKGROUND: Due to its central and strategic position in Europe and in the Mediterranean Basin, the Italian Peninsula played a pivotal role in the first peopling of the European continent and has been a crossroad of peoples and cultures since then. AIM: This study aims to gain more information on the genetic structure of modern Italian populations and to shed light on the migration/expansion events that led to their formation. SUBJECTS AND METHODS: High resolution Y-chromosome variation analysis in 817 unrelated males from 10 informative areas of Italy was performed. Haplogroup frequencies and microsatellite haplotypes were used, together with available data from the literature, to evaluate Mediterranean and European inputs and date their arrivals. RESULTS: Fifty-three distinct Y-chromosome lineages were identified. Their distribution is in general agreement with geography, southern populations being more differentiated than northern ones. CONCLUSIONS: A complex genetic structure reflecting the multifaceted peopling pattern of the Peninsula emerged: southern populations show high similarity with those from the Middle East and Southern Balkans, while those from Northern Italy are close to populations of North-Western Europe and the Northern Balkans. Interestingly, the population of Volterra, an ancient town of Etruscan origin in Tuscany, displays a unique Y-chromosomal genetic structure.

Phylogeography of human Y-chromosome haplogroup Q3-L275 from an academic/citizen science collaboration.

BACKGROUND: The Y-chromosome haplogroup Q has three major branches: Q1, Q2, and Q3. Q1 is found in both Asia and the Americas where it accounts for about 90% of indigenous Native American Y-chromosomes; Q2 is found in North and Central Asia; but little is known about the third branch, Q3, also named Q1b-L275. Here, we combined the efforts of population geneticists and genetic genealogists to use the potential of full Y-chromosome sequencing for reconstructing haplogroup Q3 phylogeography and suggest possible linkages to events in population history. RESULTS: We analyzed 47 fully sequenced Y-chromosomes and reconstructed the haplogroup Q3 phylogenetic tree in detail. Haplogroup Q3-L275, derived from the oldest known split within Eurasian/American haplogroup Q, most likely occurred in West or Central Asia in the Upper Paleolithic period. During the Mesolithic and Neolithic epochs, Q3 remained a minor component of the West Asian Y-chromosome pool and gave rise to five branches (Q3a to Q3e), which spread across West, Central and parts of South Asia. Around 3-4 millennia ago (Bronze Age), the Q3a branch underwent a rapid expansion, splitting into seven branches, some of which entered Europe. One of these branches, Q3a1, was acquired by a population ancestral to Ashkenazi Jews and grew within this population during the 1st millennium AD, reaching up to 5% in present day Ashkenazi. CONCLUSIONS: This study dataset was generated by a massive Y-chromosome genotyping effort in the genetic genealogy community, and phylogeographic patterns were revealed by a collaboration of population geneticists and genetic genealogists. This positive experience of collaboration between academic and citizen science provides a model for further joint projects. Merging data and skills of academic and citizen science promises to combine, respectively, quality and quantity, generalization and specialization, and achieve a well-balanced and careful interpretation of the paternal-side history of human populations.

Exploring the Y Chromosomal Ancestry of Modern Panamanians.

Geologically, Panama belongs to the Central American land-bridge between North and South America crossed by Homo sapiens >14 ka ago. Archaeologically, it belongs to a wider Isthmo-Colombian Area. Today, seven indigenous ethnic groups account for 12.3% of Panama's population. Five speak Chibchan languages and are characterized by low genetic diversity and a high level of differentiation. In addition, no evidence of differential structuring between maternally and paternally inherited genes has been reported in isthmian Chibchan cultural groups. Recent data have shown that 83% of the Panamanian general population harbour mitochondrial DNAs (mtDNAs) of Native American ancestry. Considering differential male/female mortality at European contact and multiple degrees of geographical and genetic isolation over the subsequent five centuries, the Y-chromosome Native American component is expected to vary across different geographic regions and communities in Panama. To address this issue, we investigated Y-chromosome variation in 408 modern males from the nine provinces of Panama and one indigenous territory (the comarca of Kuna Yala). In contrast to mtDNA data, the Y-chromosome Native American component (haplogroup Q) exceeds 50% only in three populations facing the Caribbean Sea: the comarca of Kuna Yala and Bocas del Toro province where Chibchan languages are spoken by the majority, and the province of Colón where many Kuna and people of mixed indigenous-African-and-European descent live. Elsewhere the Old World component is dominant and mostly represented by western Eurasian haplogroups, which signal the strong male genetic impact of invaders. Sub-Saharan African input accounts for 5.9% of male haplotypes. This reflects the consequences of the colonial Atlantic slave trade and more recent influxes of West Indians of African heritage. Overall, our findings reveal a local evolution of the male Native American ancestral gene pool, and a strong but geographically differentiated unidirectional sex bias in the formation of local modern Panamanian populations.

The phylogenetic and geographic structure of Y-chromosome haplogroup R1a.

R1a-M420 is one of the most widely spread Y-chromosome haplogroups; however, its substructure within Europe and Asia has remained poorly characterized. Using a panel of 16 244 male subjects from 126 populations sampled across Eurasia, we identified 2923 R1a-M420 Y-chromosomes and analyzed them to a highly granular phylogeographic resolution. Whole Y-chromosome sequence analysis of eight R1a and five R1b individuals suggests a divergence time of ∼25,000 (95% CI: 21,300-29,000) years ago and a coalescence time within R1a-M417 of ∼5800 (95% CI: 4800-6800) years. The spatial frequency distributions of R1a sub-haplogroups conclusively indicate two major groups, one found primarily in Europe and the other confined to Central and South Asia. Beyond the major European versus Asian dichotomy, we describe several younger sub-haplogroups. Based on spatial distributions and diversity patterns within the R1a-M420 clade, particularly rare basal branches detected primarily within Iran and eastern Turkey, we conclude that the initial episodes of haplogroup R1a diversification likely occurred in the vicinity of present-day Iran.

Phylogenetic applications of whole Y-chromosome sequences and the Near Eastern origin of Ashkenazi Levites.

Previous Y-chromosome studies have demonstrated that Ashkenazi Levites, members of a paternally inherited Jewish priestly caste, display a distinctive founder event within R1a, the most prevalent Y-chromosome haplogroup in Eastern Europe. Here we report the analysis of 16 whole R1 sequences and show that a set of 19 unique nucleotide substitutions defines the Ashkenazi R1a lineage. While our survey of one of these, M582, in 2,834 R1a samples reveals its absence in 922 Eastern Europeans, we show it is present in all sampled R1a Ashkenazi Levites, as well as in 33.8% of other R1a Ashkenazi Jewish males and 5.9% of 303 R1a Near Eastern males, where it shows considerably higher diversity. Moreover, the M582 lineage also occurs at low frequencies in non-Ashkenazi Jewish populations. In contrast to the previously suggested Eastern European origin for Ashkenazi Levites, the current data are indicative of a geographic source of the Levite founder lineage in the Near East and its likely presence among pre-Diaspora Hebrews.

The first peopling of South America: new evidence from Y-chromosome haplogroup Q.

Recent progress in the phylogenetic resolution of the Y-chromosome phylogeny permits the male demographic dynamics and migratory events that occurred in Central and Southern America after the initial human spread into the Americas to be investigated at the regional level. To delve further into this issue, we examined more than 400 Native American Y chromosomes (collected in the region ranging from Mexico to South America) belonging to haplogroup Q - virtually the only branch of the Y phylogeny observed in modern-day Amerindians of Central and South America - together with 27 from Mongolia and Kamchatka. Two main founding lineages, Q1a3a1a-M3 and Q1a3a1-L54(xM3), were detected along with novel sub-clades of younger age and more restricted geographic distributions. The first was also observed in Far East Asia while no Q1a3a1-L54(xM3) Y chromosome was found in Asia except the southern Siberian-specific sub-clade Q1a3a1c-L330. Our data not only confirm a southern Siberian origin of ancestral populations that gave rise to Paleo-Indians and the differentiation of both Native American Q founding lineages in Beringia, but support their concomitant arrival in Mesoamerica, where Mexico acted as recipient for the first wave of migration, followed by a rapid southward migration, along the Pacific coast, into the Andean region. Although Q1a3a1a-M3 and Q1a3a1-L54(xM3) display overlapping general distributions, they show different patterns of evolution in the Mexican plateau and the Andean area, which can be explained by local differentiations due to demographic events triggered by the introduction of agriculture and associated with the flourishing of the Great Empires.

Mitogenomes from two uncommon haplogroups mark late glacial/postglacial expansions from the near east and neolithic dispersals within Europe.

The current human mitochondrial (mtDNA) phylogeny does not equally represent all human populations but is biased in favour of representatives originally from north and central Europe. This especially affects the phylogeny of some uncommon West Eurasian haplogroups, including I and W, whose southern European and Near Eastern components are very poorly represented, suggesting that extensive hidden phylogenetic substructure remains to be uncovered. This study expanded and re-analysed the available datasets of I and W complete mtDNA genomes, reaching a comprehensive 419 mitogenomes, and searched for precise correlations between the ages and geographical distributions of their numerous newly identified subclades with events of human dispersal which contributed to the genetic formation of modern Europeans. Our results showed that haplogroups I (within N1a1b) and W originated in the Near East during the Last Glacial Maximum or pre-warming period (the period of gradual warming between the end of the LGM, ∼19 ky ago, and the beginning of the first main warming phase, ∼15 ky ago) and, like the much more common haplogroups J and T, may have been involved in Late Glacial expansions starting from the Near East. Thus our data contribute to a better definition of the Late and postglacial re-peopling of Europe, providing further evidence for the scenario that major population expansions started after the Last Glacial Maximum but before Neolithic times, but also evidencing traces of diffusion events in several I and W subclades dating to the European Neolithic and restricted to Europe.

Reconciling migration models to the Americas with the variation of North American native mitogenomes.

In this study we evaluated migration models to the Americas by using the information contained in native mitochondrial genomes (mitogenomes) from North America. Molecular and phylogeographic analyses of B2a mitogenomes, which are absent in Eskimo-Aleut and northern Na-Dene speakers, revealed that this haplogroup arose in North America ∼11-13 ka from one of the founder Paleo-Indian B2 mitogenomes. In contrast, haplogroup A2a, which is typical of Eskimo-Aleuts and Na-Dene, but also present in the easternmost Siberian groups, originated only 4-7 ka in Alaska, led to the first Paleo-Eskimo settlement of northern Canada and Greenland, and contributed to the formation of the Na-Dene gene pool. However, mitogenomes also show that Amerindians from northern North America, without any distinction between Na-Dene and non-Na-Dene, were heavily affected by an additional and distinctive Beringian genetic input. In conclusion, most mtDNA variation (along the double-continent) stems from the first wave from Beringia, which followed the Pacific coastal route. This was accompanied or followed by a second inland migratory event, marked by haplogroups X2a and C4c, which affected all Amerindian groups of Northern North America. Much later, the ancestral A2a carriers spread from Alaska, undertaking both a westward migration to Asia and an eastward expansion into the circumpolar regions of Canada. Thus, the first American founders left the greatest genetic mark but the original maternal makeup of North American Natives was subsequently reshaped by additional streams of gene flow and local population dynamics, making a three-wave view too simplistic.

Y-chromosome diversity in modern Bulgarians: new clues about their ancestry.

To better define the structure and origin of the Bulgarian paternal gene pool, we have examined the Y-chromosome variation in 808 Bulgarian males. The analysis was performed by high-resolution genotyping of biallelic markers and by analyzing the STR variation within the most informative haplogroups. We found that the Y-chromosome gene pool in modern Bulgarians is primarily represented by Western Eurasian haplogroups with ∼ 40% belonging to haplogroups E-V13 and I-M423, and 20% to R-M17. Haplogroups common in the Middle East (J and G) and in South Western Asia (R-L23*) occur at frequencies of 19% and 5%, respectively. Haplogroups C, N and Q, distinctive for Altaic and Central Asian Turkic-speaking populations, occur at the negligible frequency of only 1.5%. Principal Component analyses group Bulgarians with European populations, apart from Central Asian Turkic-speaking groups and South Western Asia populations. Within the country, the genetic variation is structured in Western, Central and Eastern Bulgaria indicating that the Balkan Mountains have been permeable to human movements. The lineage analysis provided the following interesting results: (i) R-L23* is present in Eastern Bulgaria since the post glacial period; (ii) haplogroup E-V13 has a Mesolithic age in Bulgaria from where it expanded after the arrival of farming; (iii) haplogroup J-M241 probably reflects the Neolithic westward expansion of farmers from the earliest sites along the Black Sea. On the whole, in light of the most recent historical studies, which indicate a substantial proto-Bulgarian input to the contemporary Bulgarian people, our data suggest that a common paternal ancestry between the proto-Bulgarians and the Altaic and Central Asian Turkic-speaking populations either did not exist or was negligible.

Sociocultural behavior, sex-biased admixture, and effective population sizes in Central African Pygmies and non-Pygmies.

Sociocultural phenomena, such as exogamy or phylopatry, can largely determine human sex-specific demography. In Central Africa, diverging patterns of sex-specific genetic variation have been observed between mobile hunter-gatherer Pygmies and sedentary agricultural non-Pygmies. However, their sex-specific demography remains largely unknown. Using population genetics and approximate Bayesian computation approaches, we inferred male and female effective population sizes, sex-specific migration, and admixture rates in 23 Central African Pygmy and non-Pygmy populations, genotyped for autosomal, X-linked, Y-linked, and mitochondrial markers. We found much larger effective population sizes and migration rates among non-Pygmy populations than among Pygmies, in agreement with the recent expansions and migrations of non-Pygmies and, conversely, the isolation and stationary demography of Pygmy groups. We found larger effective sizes and migration rates for males than for females for Pygmies, and vice versa for non-Pygmies. Thus, although most Pygmy populations have patrilocal customs, their sex-specific genetic patterns resemble those of matrilocal populations. In fact, our results are consistent with a lower prevalence of polygyny and patrilocality in Pygmies compared with non-Pygmies and a potential female transmission of reproductive success in Pygmies. Finally, Pygmy populations showed variable admixture levels with the non-Pygmies, with often much larger introgression from male than from female lineages. Social discrimination against Pygmies triggering complex movements of spouses in intermarriages can explain these male-biased admixture patterns in a patrilocal context. We show how gender-related sociocultural phenomena can determine highly variable sex-specific demography among populations, and how population genetic approaches contrasting chromosomal types allow inferring detailed human sex-specific demographic history.

Ancient migratory events in the Middle East: new clues from the Y-chromosome variation of modern Iranians.

Knowledge of high resolution Y-chromosome haplogroup diversification within Iran provides important geographic context regarding the spread and compartmentalization of male lineages in the Middle East and southwestern Asia. At present, the Iranian population is characterized by an extraordinary mix of different ethnic groups speaking a variety of Indo-Iranian, Semitic and Turkic languages. Despite these features, only few studies have investigated the multiethnic components of the Iranian gene pool. In this survey 938 Iranian male DNAs belonging to 15 ethnic groups from 14 Iranian provinces were analyzed for 84 Y-chromosome biallelic markers and 10 STRs. The results show an autochthonous but non-homogeneous ancient background mainly composed by J2a sub-clades with different external contributions. The phylogeography of the main haplogroups allowed identifying post-glacial and Neolithic expansions toward western Eurasia but also recent movements towards the Iranian region from western Eurasia (R1b-L23), Central Asia (Q-M25), Asia Minor (J2a-M92) and southern Mesopotamia (J1-Page08). In spite of the presence of important geographic barriers (Zagros and Alborz mountain ranges, and the Dasht-e Kavir and Dash-e Lut deserts) which may have limited gene flow, AMOVA analysis revealed that language, in addition to geography, has played an important role in shaping the nowadays Iranian gene pool. Overall, this study provides a portrait of the Y-chromosomal variation in Iran, useful for depicting a more comprehensive history of the peoples of this area as well as for reconstructing ancient migration routes. In addition, our results evidence the important role of the Iranian plateau as source and recipient of gene flow between culturally and genetically distinct populations.

Distinguishing the co-ancestries of haplogroup G Y-chromosomes in the populations of Europe and the Caucasus.

Haplogroup G, together with J2 clades, has been associated with the spread of agriculture, especially in the European context. However, interpretations based on simple haplogroup frequency clines do not recognize underlying patterns of genetic diversification. Although progress has been recently made in resolving the haplogroup G phylogeny, a comprehensive survey of the geographic distribution patterns of the significant sub-clades of this haplogroup has not been conducted yet. Here we present the haplogroup frequency distribution and STR variation of 16 informative G sub-clades by evaluating 1472 haplogroup G chromosomes belonging to 98 populations ranging from Europe to Pakistan. Although no basal G-M201* chromosomes were detected in our data set, the homeland of this haplogroup has been estimated to be somewhere nearby eastern Anatolia, Armenia or western Iran, the only areas characterized by the co-presence of deep basal branches as well as the occurrence of high sub-haplogroup diversity. The P303 SNP defines the most frequent and widespread G sub-haplogroup. However, its sub-clades have more localized distribution with the U1-defined branch largely restricted to Near/Middle Eastern and the Caucasus, whereas L497 lineages essentially occur in Europe where they likely originated. In contrast, the only U1 representative in Europe is the G-M527 lineage whose distribution pattern is consistent with regions of Greek colonization. No clinal patterns were detected suggesting that the distributions are rather indicative of isolation by distance and demographic complexities.

Bulgarians vs the other European populations: a mitochondrial DNA perspective.

To define the matrilineal relationships between Bulgarians and other European populations, we have evaluated the mitochondrial DNA (mtDNA) variation in a sample of 855 Bulgarian subjects from the mtDNA perspective. The molecular survey was performed by sequencing ∼750 bp of the control region, which resulted in 557 different haplotypes, and by a subsequent restriction fragment length polymorphism analysis to confirm haplogroup/subhaplogroup affiliation. The classification was carried out according to the most updated criteria as reported by van Oven and Kayser (Hum Mutat 30:386-394, 2009), allowing the identification of 45 mitochondrial clades. The observed pattern of mtDNA variation indicates that the Bulgarian mitochondrial pool is geographically homogeneous across the country, and that is characterized by an overall extremely high frequency of western Eurasian lineages. In the principal component analysis, Bulgarians locate in an intermediate position between Eastern European and Mediterranean populations, which is in agreement with historical events. Thus, while the Mediterranean legacy could be attributed to the Thracians, indigenous people that firstly inhabited the Balkans, the Eastern contribution is likely due to the Proto-Bulgarians originating from the Middle East and to the Slavs migrating from northeast Europe.

In search of the genetic footprints of Sumerians: a survey of Y-chromosome and mtDNA variation in the Marsh Arabs of Iraq.

BACKGROUND: For millennia, the southern part of the Mesopotamia has been a wetland region generated by the Tigris and Euphrates rivers before flowing into the Gulf. This area has been occupied by human communities since ancient times and the present-day inhabitants, the Marsh Arabs, are considered the population with the strongest link to ancient Sumerians. Popular tradition, however, considers the Marsh Arabs as a foreign group, of unknown origin, which arrived in the marshlands when the rearing of water buffalo was introduced to the region. RESULTS: To shed some light on the paternal and maternal origin of this population, Y chromosome and mitochondrial DNA (mtDNA) variation was surveyed in 143 Marsh Arabs and in a large sample of Iraqi controls. Analyses of the haplogroups and sub-haplogroups observed in the Marsh Arabs revealed a prevalent autochthonous Middle Eastern component for both male and female gene pools, with weak South-West Asian and African contributions, more evident in mtDNA. A higher male than female homogeneity is characteristic of the Marsh Arab gene pool, likely due to a strong male genetic drift determined by socio-cultural factors (patrilocality, polygamy, unequal male and female migration rates). CONCLUSIONS: Evidence of genetic stratification ascribable to the Sumerian development was provided by the Y-chromosome data where the J1-Page08 branch reveals a local expansion, almost contemporary with the Sumerian City State period that characterized Southern Mesopotamia. On the other hand, a more ancient background shared with Northern Mesopotamia is revealed by the less represented Y-chromosome lineage J1-M267*. Overall our results indicate that the introduction of water buffalo breeding and rice farming, most likely from the Indian sub-continent, only marginally affected the gene pool of autochthonous people of the region. Furthermore, a prevalent Middle Eastern ancestry of the modern population of the marshes of southern Iraq implies that if the Marsh Arabs are descendants of the ancient Sumerians, also the Sumerians were most likely autochthonous and not of Indian or South Asian ancestry.